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Two-particle correlations of direct photons were measured in central 208 Pb+ 208 Pb collisions at 158 
AGeV. The invariant interferometric radii were extracted for 100 < Kt < 300 MeV/c and compared 
to radii extracted from charged pion correlations. The yield of soft direct photons, Kt < 300 MeV/c, 
was extracted from the correlation strength and compared to theoretical calculations. 

PACS numbers: 25.75.-q,25.75.Gz 



Hanbury Brown- Twiss (HBT) interferometry provides 
a powerful tool to explore the space-time dimensions of 
the emitting source created in elementary particle or 
heavy ion collisions. Historically, such measurements 
have concentrated on pion pair correlations, but have also 
been applied to kaons, protons, and even heavy fragments 
0. Hadron correlations reflect the space-time extent of 
the emitting source at the time of freeze-out. 



The importance of direct photon Bose-Einstein inter- 
ferometry for investigation of the history of heavy ion 
collisions, especially of the very early phase,has been ex- 
tensively discussed in the literature 0, IH El 13 ■ It has 
been shown that the correlations of direct photons in dif- 
ferent ranges of the photon average transverse momenta 
reflect different stages of the collision. Therefore, photon- 
photon correlations can provide important information 
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complementary to that obtained from hadron correla- 
tions. Unfortunately, photon interferometry is faced with 
considerable difficulties compared to hadron interferom- 
etry due to the small yield of photons emitted directly 
from the hot zone in comparison to the huge background 
of photons produced by the electromagnetic decay of the 
final hadrons (primarily the neutral pion). For this rea- 
son there has been only one experimental measurement 
of photon-photon correlations in heavy ion collisions ob- 
tained at low incident energy and low photon momenta 
(Kt < 20 MeV/c) 0. In this letter we present first mea- 
surements of direct photon correlations in central ultra- 
relativistic heavy-ion collisions. 

A detailed description of the layout of the CERN ex- 
periment WA98 can be found in [7j. Here we briefly dis- 
cuss those subsystems used in the present analysis. The 
WA98 photon spectrometer, comprising the LEad-glass 
photon Detector Array (LED A), was located at a dis- 
tance of 21.5 m downstream from the 208 Pb target and 
provided partial azimuthal coverage over the rapidity in- 
terval 2.35 < y < 2.95. Further downstream, the total 
transverse energy was measured in the MIRAC calorime- 
ter. The total transverse energy measured in MIRAC 
was used for offline centrality selection. The analysis 
presented here was performed on the 10% most central 
208p b+ 208p b co ni g i ons w ith a total sample of 5.8 x 10 6 
events collected during runs in 1995 and 1996. No sig- 
nificant correlation signal was observed in the 20% most 
peripheral collision data sample of 3.9 x 10 6 events. 

In order to reject most of the hadron background, all 
showers reconstructed in the LEDA spectrometer were 
required to have a deposited energy of greater than 
750 MeV, well above the minimum ionizing peak en- 
ergy of 550 MeV. Hadron showers could be further re- 
jected by the requirement that the shower have a nar- 
row width, consistent with an electromagnetic shower in 
LEDA 7]. In addition, during the 1996 run period the 
LEDA charged particle veto was operational and pro- 
vided a shower sample of > 98% photon purity. 

The two-photon correlation function was calculated for 
each bin of the photon average transverse momentum, 
Kt = \pTi +Pt 2 1/2 as the ratio of the distribution of pho- 
ton pair invariant relative momenta, Qinv, where both 
photons were taken from the same event, to the same 
distribution but with the photons of the pair taken from 
different events. Sample correlation functions are shown 
m Fig. [U The ratio was normalized to have an equal 
number of pairs in the numerator and denominator. The 
correlation function has been fit with a Gaussian param- 
eterisation 

C 2 (Q inv ) = A ■ [1 + Aexp (-R 2 mv ■ Q} m )\ (1) 

with nomalization A, correlation strength A, and radius 
parameter Ri nv . 

There are a large number of effects which may give 
rise to small Qi nv correlations and mimic a direct photon 
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FIG. 1: The two-photon correlation function for all show- 
ers with Lmin > 20 cm and average photon momenta 100 < 
K T < 200 MeV/c (top) and 200 < K T < 300 MeV/c (bot- 
tom) fitted with Eq.0 The solid line shows the fit result in 
the fit region used (excluding the 7r° peak at Q 
and the dotted line shows the extrapolation into the low Qi nv 
region where backgrounds are large. 



pair Bose-Einstein correlation. These include: 1) single 
hadron or photon showers that are split into nearby clus- 
ters, 2) photon conversions, 3) HBT correlations from 
charged pions, or other hadrons, mis-identified as pho- 
tons, 4) residual photon correlations from ir° HBT corre- 
lations, 5) radiative decays of heavier resonances and 6) 
collective flow. 

Apparatus or analysis effects which may result in the 
splitting of a single shower into multiple clusters, or the 
merging of nearby showers into a single cluster may be 
investigated by studying the dependence of the corre- 
lation function on the relative distance L between the 
showers on the LEDA detector surface. These effects are 
expected to contribute strongly at small L and so can be 
suppressed effectively by a distance cut. Such a minimum 
cut on L introduces a lower cutoff in Q m „ 

The dependence of the correlation strength parameter 
A on the minimum distance cut I m j„ and the minimum 
invariant momentum Q m in used to define the fit region 
is shown in Fig.|2|for two Kt regions for narrow showers. 
The different symbols correspond to minimum distance 
cuts of L m i n = 20, 25, 30, and 35 cm (note that a single 
LEDA module is 4 cm in width). The results demon- 
strate that the extracted fit parameters vary strongly 
with L m i n when the low Qi nv region is included in the 
fit, a result attributed to apparatus effects and conversion 
background, but that stable results are obtained with a 
sufficiently large minimum separation distance cut, or by 
restricting the Q m „ fit region. When no charged veto 
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or narrow shape cuts are applied to the showers, stable 
results are also obtained, but with larger minimum dis- 
tance cut (or Qmin) required, consistent with the larger 
expected backgrounds. 
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FIG. 2: Comparison of A parameter fit results for differ- 
ent fit regions for 100 < K T < 200 MeV/c (top) and 
200 < K T < 300 MeV/c (bottom), calculated for narrow 
showers with different cuts on the minimum shower separa- 
tion distance: L m in= 20 cm - ■, 25 cm - A. 30 cm - T, 35 cm 
- O (same Qmin for each). 

As mentioned above, the observed correlations could 
be caused by residual correlations of charged pions, neu- 
trons, or conversion electrons misidentified as photons. 
To investigate possible contributions from non-photon 
contamination the correlation functions were constructed 
with four different identification criteria applied to the 
showers reconstructed in LEDA. These criteria have 
somewhat different photon efficiencies, which should not 
affect the photon correlation, but more importantly have 
very different levels of non-photon contamination which 
should only affect A if the contamination forms uncor- 
related background. The charged hadron contamination 
decreases from 37% and 22% to 16% and 4%, respec- 
tively for the two Kj< bins, after applying the narrow 
electromagnetic shower shape condition 0] and is negligi- 
ble after application of the charged veto condition. The 
correlation parameters extracted from these four types 
of correlation functions, corrected for contamination, are 
shown in Fig. El The consistency of the parameters ex- 
tracted with the different identification criteria indicates 
that the non-photon contribution to the observed corre- 
lation is not significant. 

The estimated systematic errors on the correlation fit 
parameters are summarized in Table [i] Besides the un- 
certainties associated with the apparatus and fit range 
discussed with respect to Fig. |2 and non-photon con- 
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FIG. 3: Comparison of parameters of correlation functions 
with different particle identification criteria: A - all clusters, 
T - narrow electromagnetic, □ - all neutral, O " narrow neu- 
tral electromagnetic (no significant result for high Kt)- 



tamination (Fig. |3J), the dependence on the fit function 
has been investigated. In addition to the Gaussian form 
of Eq. ^ an exponential form and a Student's t distri- 
bution form have also been used to fit the correlation 
functions. The Student's t distribution provides a good 
fit to the correlation for parameter n > 2. The variation 
of the fit results with n and comparison with the Gaus- 
sian fit results has been used to estimate the fit function 
error. Reasonable fit functions leading to significantly 
smaller values of A could not be found. The exponen- 
tial form gives larger A values but also exhibits a strong 
dependence of the fit parameters on fit range, indicating 
that it does not reproduce the shape of the observed cor- 
relations well. Nevertheless, if the exponential shape is 
considered with the limited fit range of Qi nv , the upper 
error bounds on A should be increased, allowing an even 
stronger direct photon signal. 

Finally, correlations might exist in the background de- 
cay photons, e.g. correlations due to collective flow, 
Bose-Einstein correlations of ir° 's, or from decays of heav- 
ier resonances. Monte Carlo simulations have been per- 
formed to estimate each of these effects. For these simu- 
lations, the transverse momentum and rapidity distribu- 
tions of the 7r°'s were taken from measurement Q. The 
effect of flow was investigated by introducing an ellip- 
tic flow pattern with a magnitude equal to that mea- 
sured for charged pions 0- Similarly, the effect of tt° 
Bose-Einstein correlations was introduced with the same 
parameters as measured for charged pions 0]. Finally, 
residual correlations due to decays of heavier resonances, 
were estimated by including all resonances having high 
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TABLE I: Summary of systematic errors on the photon- 
photon correlation parameters. 



Source 


100 <K T < 200 

A (%) R inv (%) 


200 < K T < 300 
A (%) (%) 


Apparatus 


7 


5 


16 


6 


Contamination 


17 


14 


42 


14 


Fit Function 


5 


5 


18 


6 


Fit Range 


8 


5 


26 


10 


Q inv Slope (Flow) 


2 


3 


12 


8 


Total Syst. Error (%) 


21 


17 


56 


21 


N Total Total Error ( % ) 


12 
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yield and large branching ratios for electromagnetic de- 
cay: Kg, K^, ?7, and u>. The heavy resonances were 
included based on experimental spectra where available 
and thermodynamic extrapolations otherwise. In all sim- 
ulations, the acceptance, identification cuts, and energy 
and position resolution of LEDA were applied. The sim- 
ulations indicated that elliptic flow results in the appear- 
ance of a small slope to the correlation function on the 
order of 5 • 10 -6 (MeV/c) -1 ; Bose-Einstein correlations 
of 7r°'s lead to a specific step-like correlation function, in 
agreement with analytical calculations of [5j; and resid- 
ual correlations due to decays of heavier resonances are 
found to be negligible. We have checked for flow effects 
in the data sample with highest statistics, fitting with 
a parameterization with an additional parameter for the 
slope, and found a slope parameter consistent with the 
simulations. However, limited statistics did not allow to 
extend this approach to all data, and so the final values 
of the correlation strength and radii have been corrected 
for the effect of elliptic flow. 

Averaging over the different PID criteria, we obtain 
the following correlation parameters: 

A 7 = 0.0028 ±0.0004(stat.)±0.0006(syst.) 

R{ nv = 5.4 ± 0.8(stat.) ± 0.9(syst.)fm, 

\ u = 0.0029 ±0.0007(stat.)±0.0016(syst.) 

R\l v = 5.8 ± 0.8(stat.) ± 1.2(syst.)fm 

for regions I) 100 < K T < 200 MeV/c, and II) 200 < 
K T < 300 MeV/c, respectively. 

The direct photon invariant radii can be compared to 
measurements of invariant interferometric radii of 7r~ for 
the same centrality selection and Kt region ^1 : Rinv 
7. 1 1 ± 0.22,6.91 ± 0.32, and 6.65 ± 0.30 fm at K T = 
125, 175, and 285 MeV/c, respectively. The similarity 
of the interferometric radii of direct photons and pions 
suggests that the direct photons of this Kt region are 
emitted in the late, hadron gas, stage of the collision. 

Under the assumption of a fully chaotic photon source, 
the direct photon yield N®' lrect is related to the corre- 
lation strength A and the total inclusive photon yield 
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FIG. 4: Yield of direct photons extracted from the strength 
of the two-photon correlation (closed circles) and by the sta- 
tistical subtraction method (open circles, or arrows indicating 
upper limits) Q- Total statistical plus systematical errors are 
shown. The model calculations are described in the text. 

N^ otal as 

N mrect/ N Total = 

The low pt direct photon yield has been extracted us- 
ing this expression and is presented in Fig. 0] (assuming 
Pt = (Kt))- The previously published direct photon 
yield at high transverse momenta obtained with the sub- 
traction method 7] is also shown. The measured direct 
photon results are compared with recent fireball model 
predictions ^3 ■ The calculated contributions to the total 
yield from the Quark Gluon Plasma and hadronic stages 
of the collision are shown. It is seen that the contribution 
from the hadronic gas phase dominates the direct pho- 
ton yield at small pr, with predicted yields below the 
experimental data. 

In summary, two-photon correlation functions have 
been measured for the first time in central Pb+Pb col- 
lision at 158 A GeV. The observed correlations are at- 
tributed to Bose-Einstein correlations of directly radi- 
ated photons. An invariant radius of about 6 fm is ex- 
tracted, comparable to that extracted for pions of simi- 
lar momenta, and the correlation strength parameter was 
used to extract the yield of direct photons at pt < 300 
MeV/c. The extracted yield exceeds theoretical expecta- 
tions which attribute the dominant contribution in this 
Pt region to the hadronic phase. 
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